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SECTION 1 
INTRODUCTION 


This manual describes the computer program written 
in support of the problem to determine Aerodynamic 
Influence Coefficients (AIC's) on parallel interfering wings 
(See Ref. 1. ) 

The information presented here is geared to the 

programmer. It is sufficient to fully describe the 

* 

program logic and the required peripheral storage. 

Figure 1 gives an overview of the entire program, and 
is the basis for the control program. All element gen- 
erated information is stored externally to reduce core 
storage. A separate section is devoted to the develop- 
ment of these external files so that I/O time may be 
optimized through efficient buffer description. In- 
dividual subroutine write-ups are presented along with 
the complete Fortran source listing. Program alterations 
(data storage) are also discussed. 
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Figure 1.1 Computer Program Overview 
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Figure 1.1 (Continued) 
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SECTION 2 

EXTERNAL FILE STRUCTURE 

This program utilizes eight (8) files during the course of 
operation. The delivery version of the AIC / INT program comes 
with the following variable names and real unit designations, any 
or all of which the programmer may alter: 


Unit 

Name 

Unit ID 

Type 

Usage 

NTAPE 


1 

(S,P) 

Velocity Potentials 

LTAPE 


2 

(S,P) 

Pressure Arrays 

MTAPE 


3 

(S,P) 

Downwash Potentials 

KTAPE 


4 

(S) 

Scratch 

SYSIN 


5 

(S) 

Standard Card Input ( 8 0 Col Card) 

SYSOUT 

, OUT 

6 

(S) 

Standard Line Printer (132 Char/Line) 

MREST 


8 

(S) 

Restart Only, Store MTAPE 

TYSIN, 

TYSOUT 

9 

(S) 

Store Card Images from SYSIN; 
also Store Title Cards 

where 

S = Scratch , P = 

Permanent . 



If the restart feature is desired then units so designated (P) must 
be permanent files. 

(NOTE: As units 5 and 6 are standard input and output, they will 

not be considered in the following discussion.) 

Unit TYSIN=TYSOUT has a logical record length (LRECL) of 20 words 
(@ 4-literal/word ) . 
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Unit NTAPE uses LRECL = 14 words only. 

Unit LTAPE, KTAPE uses LRECL = variable 

Unit MTAPE uses both o£ above. 

All complete, complex, matrices are stored on tape by Subroutine 
TINOUT. These arrays are stored column-wise. Each is preceded by 
a header, whose first word is "-10" and a trailer, whose first word 
is "-20". Element type information (LRECL = 14) has only the end 
of information defined by a trailer, whose first word is "-20". 

Figure 2 J. gives an overview of the information stored on specific 
units . 
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Those areas marked ** indicate that an array was written. by 
Subroutine TINOUT. 






SECTION 3 

S UBROUTINE WRITE-UPS 

The AIC / INT program encompasses 25 subprograms, each with 
a unique sequence number. Columns 73 , 7^, 75 contain the 
"DECK" number and columns 76 thru 80 contain the card sequence 

number for that subprogram. The first card is always 00010 

with successive increments of 10. 

Included with each description is a statement declaring the 
size of the program. This number is intended as a guide only, 
as it reflects the storage requirement of an IBM/ 360 / 65 , FORTRAN IV 
G level 19 compiler. 

In this manual the subroutine writeups are presented in 
alphabetical order. These are cross-referenced by deck number 
on the following page. 
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SUBROUTINE WRITE-UPS 


DECK NO. 


014 
012 
013 
000 
021 
022 
020 

024 
023 

003 

005 
016 
008 

025 
002 

009 
001 

006 

010 

015 
Oil 

004 
018 
019 
017 
000 
007 


NAME 

ASSDWN 

ASSPOT 

ASSY . 

BASIC 

CMATPR 

CMINV 

CMPRD 

CORED 

CRASH 

DATA . 

DIAPH3 

DISPLA 

DWASHY 

ENDJOB 

INPUT 

FAKTOR 

MAIN . 

MASTER 

ORIENT 

OUTPUT 

PSI . . 

SIZE . 

STORE 

TAPES 

TINOUT 

TRINT 

TRINTX 



1. Subroutine Name : ASSDWN 

2 . Purpose : 

Assemble downwash contribution to the potentials, condense and 
compute the final velocity potential distribution across the 
wing surface. 

3. Equations and Procedures : 

Downwash influence coefficients at nodes of wing 1 due to the 
presence of wing 2 are written as 

NJ ' [v. w j V.J 

where the subscripts w and d denote the source strength situated 
on the wing and diaphragm respectively. Similarly, for wing 2 
due to wing 1, 



The solution for the source strengths on the upper surface of 
wing 1 (GJ U ) and the lower surface of wing 2 (a; I ) is. computed 
from 
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The array P is computed in subroutine ASSPOT and P2/Pt 
is the input downwash computed in subroutine master. 


The total source strengths across the wing are given by: 



Input Arguments : 


INDWN 

NP1 

NP2 

ND1 

ND2 

NP 


input downwash array 
total no. nodes (wing only), wing 1 
total no. nodes (wing only), wing 2 
total no. nodes (diaphragm only) wing 1 
total no. nodes (diaphragm only) wing 2 
total no. nodes = NP1 + NP2 
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Input Arguments : (cont'd) 

ND = total no. nodes = ND1 + ND2 

NP2ND2 = total no. nodes - wing 2 = NP2 + ND2 

NP1ND1 = total no. nodes wing 2 = NP1 + ND1 

NMOD = no. mode shapes 

MTAPE = unit defining downwash potentials 

KTAPE = unit defining velocity potentials 

Output Arguments : 

All input, unmodified 

Error Returns : None 

Calling Sequence : 

CALL ASSDWN (W12, W21, A, WT, WORK, T, INDWN, NP, NP1 , NP2 , 

ND1, ND2, ND, MTAPE, KTAPE, NMOD, NP2ND2 , NP1ND1 
C) 

Input Tapes : 

MTAPE = down wash potentials 

KTAPE = condensed velocity potentials 

Output Tapes : None 

Scratch Tapes : None 

Storage Required : 

547^ bytes (1369 words) 

Subroutine User : MAIN 

Subroutine Required : 

CMATPR, TINOUT, STORE, CMPRD, CMINV 

Remarks : 

W12, W21, A, WT, WORK, T. C 

are all complex arrays, used as work space in subroutine. 



1. Subroutine Name: 


ASSPOT 


2 . Purpose : 


Assemble from element information the condensed velocity 
potential array. 


3 . Equations and Procedures : 




I, 


condensed 


_ ^WW _ [ ^WD 


T 


4 . Input Arguments : 


PHI 

PHIXX 

PHI 2 2 

NP1 

NP2 

ND1 

ND2 

NTAPE 

KTAPE 

MTAPE 


work space 
work space 
work space 

Nbr grid points, wing 1 
Nbr grid points, wing 2 
Nbr grid points diaphragm 1 
Nbr grid points diaphragm 2 
tape logical unit nbr 

it ii ti ii 

ii it ii it 


5 • Error Returns : None 


6 . Calling Sequence : 

Call ASSPOT ( PHI, PHIXX, PHI22 ,T,NP1 ,NP2 ,ND1 ,ND2 , NTAPE, KTAPE, MTAPE) 

7 . Input Tapes : UNIT = NTAPE 

8. Output Tapes : UNIT = KTAPE 

UNIT = MTAPE 

9 • Storage Required : 

2750 bytes (688 words) 
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ASSPOT 
Page 2 


10 . subroutine User : 
main 

11. s ubroutine Required : 

ASSY, CORED, TINPUT, CMINV, CMATPR, CRASH, CMPRD 


12. Remarks : 
None 




1. Subroutine Name: 


ASSY 


2. Purpose : 

Assemble a large matrix from individual element information 
stored on tape . 


3. Equations and Procedures : 


Records on input are exactly 14 words long. 


Word (1) 
Word (2) 
Word (3) 
Word (4) 


row position 

col position for word ( 9) 
col position for word (11) 
col position for word ( 13 ) 


4 . Input Arguments : 

PHI = Array to be assembled, previously cleared 

N = Row dimension of PHI 

NP1 = Nbr grid points, wing 1 

NP2 = Nbr grid points, wing 2 

ND1 = Nbr diaphragm points, wing .1 

NTAPE = Unit (tape) 

NSET = Assembly branch (1, 2, 3, or 4) 


5. Output Arguments : 

PHI = Assembled array. Because of program logic and/or 
problem influence, array may be partially empty. 


6. Error Returns : 

End-Of-File condition stops program with a value "7" 

7 . Calling Sequence : 

Call ASSY (PHI, N,NP1,NP2,ND1, NTAPE, NSET) 

8. Input Tapes : UNIT = NTAPE 

9. Storage Required : 2822 bytes (705 words) 

10. Subroutine User : ASSPOT 

11. Subroutine Required: None 
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NAME COMMON: /BASIC/ 

This common statement appears in all subroutines. It 
contains all the basic control data required for operation and 
control of the total program. 

The total length of this common is 100 words. Words 1 
through 15 apply only to Wing 1; words 16 through 30 apply only 
to Wing 2. The array is defined as follows: 


Word 

/Variable 

Meaning 

1,16 

NEL 

Total # elements on planform 

2,17 

CAPPA 

Frequency (pure imaginary - stored as real) 

3,18 

MLlNE 

# grid. points on mach line 

4,19 

NMOD 

# mode shapes 

5,20 

NP 

# grid points (wing only) 

6,21 

NELW 

# elements (wing only) 

7,22 

NGP 

Total # grid points on planform 

8,23 

ND 

# grid points diaphragm 

9 , 24 

NFLD 

# elements in diaphragm 

10,25 

SYMC 

Symmetry factor 

11,26 

NDIV 

# divisions diaphragm 

12,27 

REF 

Reference length 

13,28 

NLND 

# leading edge grid points 

14,29 

THETA 

Wing angle 

15,30 

- 

N/A 

30-45 


N/A 

46 

PI 

PI = 3.1415 etc 

47 

NP 

NP1 + NP2 . 

48 

- 

NP1 + ND1 

49 


NP2 + ND2 

50 

ND 

ND1 + ND2 

51 

- 

NP * NP 

52 

- 

ND * ND 

53 

NPND 

NP * ND 

54 

KTAPE 

Tape, Unit =4 

55 

LTAPE 

Tape, Unit = 2 

56 

MTAPE 

Tape, Unit =3 
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Word 

Variable 

Meaning 

57 

NTAPE 

Tape , Unit = 1 

58 

SYSOUT 

Line printer 6 

59 

SYSIN 

Card reader 5 

60 

TYSIN 

Temporary card reader 9 

61 

KPRINT 

MAIN 

62 

KPRINT 

INPUT, DATA, DIAPH3 

63 

KPRINT 

MASTER, ORIENT 

64 

KPRINT 

TRINTX 

65 

KPRINT 

PSI 

66 

KPRINT 

DWASHY 

67 

KPRINT 

ASSPOT, ASSY 

68 

KPRINT 

ASSDWN 

69 

KPRINT 

CMINV 

70 

KPRINT 

TAPES (22) 

71 

KPRINT 

STORE, TINOUT 

72 

KPRINT 

CM PRD 

73 

- 

N/A 

74 


N/A 

75 

IWING 

Master, receiving element is on 

76 

KWING 

Master, influencing element is on 

77 

L 

Master, receiving element 

78. 

M 

Master, influencing 

79 

IN 

Master, node point # 

80 

SWITCH 

Master, branch control 

81 

ISW 

Master, branch control 

82 


N/A 

83 

- 

Punch option (yes/0) 

84 

IPUN 

Punch Unit =7 

85 

EPS 

eps . 

86 

RESTART 

NO = 0, YES = -1 
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Word 


Variable 


Meaning 


87 

MACH 

Mach number 

88 

VEETA 

3 = V MACH 2 -1 

89 

NTITLE 

# title cards 

'90 

ISOLAT 

Extrapolation; Yes =1, No = 0 

91 

NEWPTS 

# new grid points for extrapolation 

92 

MU 

Mach angle 

93 

D 

Stagger distance 

94 

- 

N/A 

95 

- 

N/A 

96 

- 

N/A 

97 ! 

DEPS 

Convergence criteria, PSI (sum) 

98 

- . 

N/A 

99 

NERR 

Errors generated 

100 

RUN 

Yes = 1, No =0 
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CMATPR 


1 . Subroutine Name: 


2. Purpose : 

To display a complex array of order N,M in a formal fashion. 


3. Equations and Procedures : 

Total number of columns to be printed is broken into integral 
groups of 4 and displayed as rows 1 to N, columns I to J, . 
when (J-I + l)-4 


4. Input Arguments : 

A = Complex array, order N#M exactly 
N = Number of rows in "A" 

M = Number of columns in "A" 

NAME = Literal data, 4 characters long 

5- Output Arguments : None 

6. Error Returns : None 

7. Calling Sequence : Call CMATPR ( A,N,M, 4HNAME) 

8. Input Tapes : None 

9- Output Tapes : None 

10. Scratch Tapes : None 

11. Storage Required : 1196 bytes (299 words) 

12. Subroutine User : DISPLA, OUTPUT, ASSDWN, ASSPOT 

13* Subroutine Required : None 

14 . Remarks : 

The literal data stored in NAME is a control for header infor- 
mation to be displayed from this subroutine. If NAME = 4HN0NE, 
no header information is displayed. Otherwise, NAME ,N,M are 
displayed at top of the array. 
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1. Subroutine Name : CMINV 

2. Purpose : Invert a matrix 

3 . Equations and Procedures : 

The standard Gauss-Jordan method is used. The determinant 
is also calculated. A determinant of zero indicates that 
the matrix is singular. 

4 . Input Arguments : 

A = input matrix; complex 
N = order of "A" 

5- Output Arguments : 

A = inverse of input matrix ; complex 

6 . Error Returns : None 

7 . Calling Sequence : 

CALL CMINV (A, N) 

8 . Storage Required : 3484 bytes (871 words) 

9. Subroutine User : ASSDWN, ASSPOT 

10. Subroutine Required : None 

1 1 . Remarks : 

a. Input array must be stored as a general matrix 

b. A singular array terminates program with a message. 





1 . 


Subroutine Name: 


CM PRD 


2 . Purpose : 

Special purpose complex matrix multiplication. 

R(N,L) = A(N,M) * B(M,L) 

3- Equations and Procedures : 

m 

r(i,J) = a ( i , k ) * b(k,J) • 

All arrays are treated as singularly subscripted, stored 
column wise and fully packed. This is a row into column 
product . 

A . Input Arguments : 

A = pre-multiplier array, order (N,M): complex 
B = post-multiplier array, order ( M , L ) : complex 
N = row dimension of "A" 

M = column dimension of "A" 

L - column dimension of "B" 

5. Output Arguments : 

R = resultant array, order (N,L) : complex 

6. Error Returns : None 

7 . Calling Sequence: 

CALL CMATPR (A,B,R,N,M,L) 

8. Storage Required : 728 bytes (182 words) 

9 • Subroutine User : 

OUTPUT 
ASSDWN 
ASSPOT ' 

10. Subroutine Required : None 

11. Remarks : All three arrays must occupy separate locations 
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1. Subroutine Name : CORED 

2. Purpose : 

Purpose is to force a CORE DUMP of main storage as an aid to 
locating the reason for calling this routine. 


! 

3. Equations and Procedures : Self-explanatory 

4 . Input Arguments : None 


5 . 

Output Arguments: 

6 . 

Error Returns: 

7 . 

Calling Sequence: 

8 . 

Storage Required: 

9 . 

Subroutine User: 

10 . 

Subroutine Required: 


None 

None 

CALL CORED 
308 bytes (77 words) 
CMINV, ASSPOT, DIAPH3 
None 


11. Remarks : 

**NOTE** This subroutine must be modified to fit the machine 
on which it is used. The minimum requirement of this sub- 
routine is to CALL EXIT. 
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1 . 

2 . 

3- 

4 . 

5. 

6 . 
7- 
8 . 


Subroutine Name : CRASH 

Purpose : 

To display a full page message, "ERROR HALT", indicating 
that job has terminated abnormally. 

Equations and Procedures : Self-explanatory 

Callin g Sequence : CALL CRASH 

Storage Required : 884 bytes (221 words) 

Subroutine User : CMINV, ASSPOT, DIAPH3 , DATA, INPUT 

Subroutine Required : None 

Remarks: None 
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1. Subroutine Name: 


DATA 


2 . Purpose : 

Read, sort and store all input data. Perform tests where 
possible and track input errors. Display input data. 

3. Equations and Procedures : Self-explanatory 

4 . Input Arguments : 

IW = kode (1 or 2) to identify call number 

5 . Output Arguments : 

XYZ = grid point coordinate data 

IBLN = element Boolean 

XMOD = mode shape data 

LND - leading edge grid points 

TITLE = work space for storing title cards 

6. Error Returns : KONTRL(99) > 0 means stop 

7. Calling Sequence : 

Call DATA (XYZ, IBLN, XMOD, LND, TITLE, IW) 


8. 

Input Tapes: TYSIN 

(temporary card image 

9. 

Output Tapes: None 


10. 

Scratch Tapes : None 


11. 

Storage Required: 

3462 bytes (866 words) 

• 

C\J 

i — 1 

Subroutine User: 

INPUT 

13. 

Subroutine Required: 

CRASH, EXIT 

14. 

Remarks: None 
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1. Subroutine Name : DIAPH3 


2. Purpose : 

Generate grid point coordinates and element Boolean for the 
diaphragm region. 


3 . Equations and Procedures : 

Each new grid point is 
generated from a pair of 
corresponding points 
along the leading edge. 
Since the angle 0 and u 
are known, the length IJ 
can be computed. It is 
then a simple matter to 
compute the X,Y coordinates 
of J . 


4 . Input Arguments : 



XY 

IBLN = 

LND 

IAR 

NLND = 

NEL 

NGP 

NDIV = 
MACH = 
MU 
OUT 
PI 


X,Y coordinates array 
existing Boolean (element) 
leading edge grid point numbers 
work space 
number of LND 

number of elements on wing 

number of grid points on 

number of divisions along mach line 

mach number 

mach angle 

SYSOUT value 

PI 


5. Output Arguments : 


LDNW 

NLDNW 

NEL 

NGP 

NERR 


mach line grid points (excluding lead point) 

number of LDNW 

total number of elements 

total number of grid points 

error condition 


6. Error Returns: 


NERR > 0 implies negative grid points generated; either logic 
has failed or input data in error. 
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Calling Sequence : 

Call DIAPH3 (XY,IBLN,LND,LDNW,IAR,NLND,NLDNW,NEL,NGP 
NELD1 , NDIV, MACH , MU , OUT , PI ,NERR ) 

Storage Required : 

3^62 bytes (866 words) 

Subroutine User : 

INPUT 

Subroutine Required : 

CRASH, CORED 

Remarks : 


None 


1 . 


Subroutine Name: 


DISPLA 


2. Purpose : 

Display a portion of fixed length complex array. 

3. Equations and Procedures : 

That portion of the array to be displayed is transferred to 
a. fixed length array and then the standard complex matrix 
printout is called. 


4. 


Input Arguments : 


INDWN 

WORK 

NP1 

NP2 

NP 

NMOD 


=■ complex array, (100, 6), to be displayed 
= complex array, (NP,NM0D) work space 
= not used 
= not used 

- NP1 + NP2 = number of rows of INDWN to be displayed 
= nbr. of columns of INDWN to be displayed 


5- Output Arguments : None 

6. Error Returns : None 

7. Calling Sequence : 

CALL DISPLA . ( INDWN , WORK , NP1 , NP2 , NP , NMOD ) 

8. Storage Required : 690 bytes (173 words) 

9 . Subroutine User : MAIN 

10. Subroutine Required : CMATPR 

11. Remarks: None 




1. Subroutine Name: 


DWASHY 


2. Purpose : To evaluate the following integral across a 

triangular element using Gaussian quadrature 
technique . 


1 - 


r 


- e -4 (*■-*■) y 

j- - ^ V " ?*■) 'v \r ' \ 

* Y(^*) 

CXr - fj ' 


el 


)0- 


3. Equations and Procedures : 

The equations and procedures required to compute the above 
integral are the same as those used in Subroutine TRINTX. 

The only difference here is that this is a single integral. 

4 . Input Arguments : 

Yl, Y2, Y3 = ordered set of Y-coord values of the triangle 

(smallest to largest) 

Z = Z-coord of triangle (Zl+Z2+Z3)/3 

Al, A2, A3 = slopes of sides of triangle 


Bl, B2, B3 = Y-intercept of slopes 

Case = (1 or 2) relative location of the base of 

triangle 

5. Output Arguments : None 

6. Error Returns : None 

7. Calling Sequenc e: 

CALL DWASHY (Yl, Y2, Y3, A12, A13, A23, B12, -B13, B23, CASE, 
Z) . 

8. Storage Required : . 2970 bytes (7^3 words) 

9 • Subroutine User : MASTER 
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10. Subroutine Req ui red : 


PSI 


11. Remarks: 




A typical element Is shown here to indicate what points 
are chosen for integration. 
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1 . 


Subroutine Name: 


ENDJOB 


2 . Purpose : 

Display a full page title stating that program has termi- 
nated normally. 


3. Equations and Procedures : Self-explanatory 


4 . Input Arguments' : 


None 


5 . Output Arguments 


None 


6. Error Returns: 


None 


7 . Calling Sequence 


Call ENDJOB 


8 . Storage Required : 


820 bytes (205 words) 


9 • Subroutine User : 


MAIN 


10. Subroutine Required : 


None 


11. Remarks: 


Page centering is accomplished by adjustment of the prefix 
N in line 1 of the format -statement 1132, 
e.g. (... N ( 1H / ) . . . ) 
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1 . 


Subroutine Name 


INPUT 


2. Purpose : 

Read input stream and store on temporary data set. Display 
card images as read. Control reading and storing (Sub. DATA) 
of data and diaphragm generation (Sub.DIAPH3) . 


3- Equations and Procedures : Self-explanatory 

4. Input Arguments : None 


5 . Output Arguments : 


XYZ 

IBLN 

XMOD 

LND 

MLINE 

TITLE 


grid point coordinate data 
element Boolean 
mode shape data 

leading edge grid point numbers 
mach line grid point numbers 
work space 


6. Error Returns: None 


7- Calling Sequence : 

Call INPUT (XYZ, IBLN, XMOD, LND, MLINE, TITLE) 


8. Scratch Tapes : TYSIN = TYSOUT 


9. Storage Required: 4020 bytes (1005 words) 


10. Subroutine User: MAIN 


1 1 . Subroutine Required : 

CRASH 

DATA 

DIAPH3 


12. Remarks: None 
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1. Subroutine Name: 


FAKTOR 


2. Purpose : To compute the total number of elements having a 

common node point (wing only) and return this value as 
a weighing factor for each node point. 

3* Equations and Procedures : 

Each element contains three node points stored in an sequential 
array. It thus ' requires only the comparison of a given node 
point number, say i, to all node points in the list. The 
number of appearances is counted and the reciprocal of this 
value is the weight factor. 

4 . Input Arguments : 

IBLN = Boolean array of node points vrs element no. 

NEL = Total number of elements (wing only) 

NGP = Total number of grid points (wing only) 

5. Output Arguments : 

FAC = Weight factor per node point 

6. Error Returns : None 

7. Calling Sequence : 

CALL FAKTOR (FAC , IBLN, NFL, NGP ) 

8. Input Tapes : None 

9* • Output Tapes : None 

10. Scratch Tapes : None 

11. Storage Required : 

552 bytes (138 words) 

12. Subroutine User : MASTER 

13. Subroutine Required : None 

14. Remarks: None 
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1 . Subroutine Name : 


MASTER 


2. Purpose : This routine is the master logic control program 

for generation of element to element interaction data 

3. Equations and Procedures : 

The logical flow of this routine is shown in the accompanying 
flow chart. In addition to this flow chart, pressure matrices 
and' input downwash matrices are computed as follows: 

Pressure matrices for the basic elements are: 


[A] -■ [£] + cs'-o 

[ft] - 



2 A - Aei 


^ I ^ V 


y, 

1 


» 1 


-n. x -- ()w, ),(y.->0] 



Input downwash, at each node 

(U) . jn, . |X 
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4. Input Arguments : 


XMOD = mode shape data 

XYZ = coordinate data 

MLINE = mach line node points 

IBLN = element . boolean array 

LND = leading edge node points . 

5 • Output Arguments : 

INDWN = input downwash array 

6. Error Returns : None 

7 . Calling Sequence : 

' I 

CALL MASTER (INDWN, XMOD, XYZ, MLINE, IBLN, PM, PTUSED, LND) 

8. Input Tapes : None 

9. Output Tapes : 

NTAPE = (Unit=l) , velocity potentials 

MTAFE = (Unit=3) , downwash potentials 

LTAPE = (Unit=2) , pressure matrices 

10. Scratch Tapes : 

11 • Stor a ge Required : 

7428 bytes, (1857 words) 

12. Subroutine User : MAIN 

13- Subr outine Required : 

FAKTOR , DWASHY., TRINTX 
ORIENT , TRINOUT 

14. Remarks : None 
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Subroutine Name 


ORIENT 


Purpose : 

To determine the orientation of a 2-dimensional triangle in 
space and compute the slopes and intercepts of all 3 sides. 

Equations and Procedures : 

a) arbitrarily choose an order for starting 

b) through a logical comparison determine the new order 

c) so ordered, the slopes and Intercepts are computed by 
y = mx + b 

Input Arguments: 


X = array of 3 -x coordinates 
Y = array of 3 -y coordinates 

Output Arguments : 


XI, X2, X3 
Yl, Y2 , Y3 
A12,A13,A23 
B12,B13,B23 
CASE 


ordered x-values, low to high 
corresponding y values 
slopes of lines 
intercepts of lines 

normally 1 but set = 2 If base of triangle 
is parallel to y-axis 


Error Returns : None 


Calling Sequence : 

CALL ORIENT ( X, Y, A12 , A13, B12, B13, CASE, 

XI, X2 , X3 j A 2 3 j B23 , Yl, Y2 , Y3) 

Subroutine User : Master 

Subroutine Required : None 

Remarks : None 



1. Subroutine Name ; OUTPUT 

2 . Purpose : 

Display the major results of the total program. Some final 
calculations are performed prior to display. 

3 • Equations and Procedures : 

[PRESSURES] = [PRESSURES INPUT] * [VELOCITY POTENTIALS] 
[CEN. FORCES] = [MODESHAPES] # [PRESSURES] 


4 . 


Input Arguments : 


PHI 

PM 

XMOD 

PRES 

GENF 

NP 

NMOD 

LTAPE 

HEY 


velocity potentials (NP # NMOD) 

work space (pressure matrix storage) 

mode shapes (50 * 6 * 2) 

work space (final pressures) 

work space (generalized forces) 

order of wing for display 

Nbr. mode shapes used 

Unit definition for pressures 

(1 or 2) describes wing for display 


5- Output Arguments : None 

6. Error Returns: None 


7. 

8 . 


9- 

10 . 

11 . 

12 . 


13. 


14. 


Calling Sequence : 

CALL OUTPUT (PHI, PM, XMOD, PRES, GENF,NP, NMOD, LTAPE, KEY) 
Input Tapes : 


UNIT = LTAPE 
UNIT = TYSIN 


Output Tapes : 

Scratch Tapes : 
Storage Required : 

Subroutine User: 

Subroutine Required : 

Remarks : 


None 

None 

2090 bytes (525 words) 
MAIN 

TINOUT, CMPRD, CMATPR 
COMMQN/BASIC/KON ( 100 ) 
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1. Subroutine Name : PSI 

2. Purpose : Evaluate the function described below 

3. Equations and Procedures :- 

Yu*) • 5>"(.”v » 

where the J n are Bessel functions of integer order. 
Sufficient terms are taken to insure error in the sum 
-0.0005. ■ 

4 . Input Arguments : 

XJ = x-coordinate 

YJ = Y-coordinate 

g - g-coordinate 

5. Output Arguments : 

PSI = computed value "Y(u,v) 

6 . Error Returns : 

IER Meaning 

1 max terms in sum (N) was >40 

2 max terms in sum (N) was > 100 

3 argument of Bessel function was negative 

4 convergence criteria on sum was not met 

All messages non-fatal. Final results should be suspect 

7. Calling Sequence : 

Call PSI (XJ,YJ,g,SEI) 

8 . Storage Required : 

3736 bytes (944 words) 
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1 . 


Subroutine Name: SIZE 


2. Purpose: To compute a) the starting locations of all variable 

length arrays used during analysis, b) test for 
maximum work space exceeding that alloted to program. 
If b is true, message is returned to, calling program. 

3. Equations and Procedures : 

Using input data concerning the current problem, all work space 
for each Subroutine to be called is computed by successive 
additions of each array length to a grand total. Any array 
being shared by more than one. subroutine is located near the top 
of the work space. When all space computed, grand totals are 
compared to the alloted maximums. 


4. 

Input Arguments: 

KON = 

COMMON/BASIC/KONTRL , 



IS 

contains maximum array sizes and 
maximum lengths of work space 



KEY = 

1 or 2 for branching 

5. 

Output Arguments: 

IS 

contains starting locations of all 




work space arrays 


6. Error Returns : If maximum size(s) exceeded, KON(IOO) = 0 is 

returned as an. error flag. 

7. Calling Sequence : CALL SIZE (KON, IS, KEY) 

8. Storage Required : 1598 bytes ( 4 0 0 words) 

9. Subroutine User : MAIN 

10. Subroutine Required: None 
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Subroutine Name ; STORE 
Purpose ; 

To insert a submatrix into a master array 
Equations and Procedures : 

Program is capable of inserting elements as follows: 


(a 

= master 

array ) 

a. 

a ( i , j ) 

= b(k,l) 

b . 

a ( i , j ) 

= -b (k,l ) 

c . 

a ( i , j ) 

= a ( i , j ) - b(k,l) 


Each method is selected by a control variable "KODE" 
(see input argument list). 


Input Arguments: 


WT 

= 

master array, complex 



W 

= 

submatrix, complex 



NP . 

= 

row dimension of "WT" 



NN 

= 

row dimension of "W" 



Ml 

s= 

starting subscript (rows) 

in 

"WT" 

LI 

= 

starting subscript (cols) 

in 

it yip" 

M2 

= 

stopping subscript (rows) 

in 

"WT" 

L2 

= 

stopping subscript (cols) 

in 

»»yip» 

KODE 

1 “ 

select mode of storage 




1 = type a above 

2 = type b above 

3 = type c above 


Output Arguments : None 

Error Returns : None 

Calling Sequence : 

Call STORE (WT, W y NP , NN, Ml, LI, M2 , L 2, KODE) 

Storage Required : 1390 bytes (347 words) 

Subroutine User : ASSDWN 

Subroutine Required : None 

Remarks : 

It is assumed that before KODE = 3 is used that "WT" array 
contains valid data. 



1 . 

Subroutine Name: 

TAPES 


2. 

Purpose : 




Special purpose subroutine to display 3 (threee) tapes used 
in the AIC/INT program. These tapes store all element related 
information so generated, and are the basis of the RESTART 
phase of the program. 

3. 

Equations and Procedures: 

Self-explanatory 


4. 

Input Arguments: 




C = array of complex storage, 2500 words will handle 
maximum case generated. 

the 

5. 

Output Arguments: 

None 


6. 

Error Returns: 

None 


7- 

Calling Sequence: 
CALL TAPES (C) 



8. 

Input Tapes : 




UNIT =1 
UNIT =3 
UNIT =2 



9. 

Output Tapes: 

None 


10. 

Scratch Tapes: 

None 


11. 

Storage Required: 

984 bytes (246 words) 


12. 

Subroutine User: 

MAIN 


13- 

Subroutine Required: 

C MAT PR, TINOUT 


14. 

Remarks : 




This subroutine was designed specifically for "DEBUG" purposes. 

#*WARNING*# placement of call to this subroutine is critical, 
as all tapes are left "REWOUND". This is an abnormal situation 
at any point in program except the end. 
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1 . 


TINOUT 


2 . 

3. 


4. 


5. 

6 . 

7. 

8 . 
9- 

10 . 

11 . 

12 . 


Subroutine Name: 


Purpose : 

Read/write a complex array of order (L,M) from/to tape. 


Equations and Procedures : 


Each array Is placed 

a. Identification 

b. Array follows, 

c. Identification 


on tape as follows: 
header is written first; 
stored by columns; 
trailer is written last; 


HEADER : 

Word 1 2 3 4 



7 


8 thru 14 


-10 


L M N A M E 


0 


TRAILER: Word 1 = -20, 

Remainder record the same as header 


Input Arguments : 


X 

L 

M 

NAME 

10 

KTAPE 


array, complex, to be written to tape 
row dimension of X 
column dimension of X 

4 literal characters describing the name or- array X 
control for action taken 

unit where information will be read/written 


Output Arguments : 

X = array, complex, read from tape 


Error Returns : 

Calling Sequence : 
CALL TINOUT (X, L, M 

Input Tapes : 

Output Tapes : 

Storage Required : 

Subroutine User: 

Subroutine Required: 


None 

NAME, 10, KTAPE) 

UNIT = KTAPE (10=2) 

UNIT = KTAPE (10=1) 

1722 bytes (431 words) 
MASTER, ASSP0T, ASSDWN, 
CORED 


OUTPUT 
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Remarks : 


Four error conditions may arise because of failure in other 
parts of the program (AIC/INT). These errors will cause 
termination of the program thru the call to CORED. All 
four errors involve failure to match input information 
(thru argument) with that read from the tape. 

ERROR NBR MEANING 

First record not a header 

Name of array requested does not match that 
found in header 

Order of array requested disagrees with that 
found in header 

Last record read was not a trailer 


1 

2 

3 

4 


Any of these messages indicate the program logic has failed. 



NAME COMMON: /TRINT/ 


The total length of this common block is 56 words. It appears In 


only four 
variables 
follows : 

subroutines : 
in the block 

MASTER, 
(in order 

TRINTX, DWASHY, and PSI. The 
of appearance) are defined as 

Variable 

Type 

Dimension 

Definition 

FPOT 

C 


3 

Velocity potentials 

DWN 

C 


3 

Downwash potentials, X-integration 

DWNI 

C 


3 

Downwash potentials, Y-integration 

UU 

R 


5 

Pivot points, Gaussian integration 

CC 

R 


5 

Weight Factors, Gaussian integration 

1 

XI ' 

R 


3 

X-coord, influencing element 

YI 

R 


3 

Y-coord, influencing element 

ZI ' 

R 


3, 

Z-coord, influencing element 

1 ■ 

T 

R 


9 

Transformation, influencing element 

XR 

R 


1 

X-coord, receiving element 

YR 

R 


1 

Y-coord, receiving element 

ZR 

R 


1 

Z-coord, receiving element 

XK 

R 


1 

Frequency 

ZMACH 

R 


1 

Mach number 

BEETA 

R 


1 

beta 

PIE 

R 


1 

Pi 

XKOM 

R 


1 

Frequency / mach no. 

USQRD 

R 


1 

beta * beta 

C10P.IV 

R 


1 

1 / (pi * beta) 



1. Subroutine Name: 


TRINTX 


2. 



Purpose: To evaluate the following integrals across a 

triangular element using Gaussian quadrature 
technique. 



SX&r-V - nft.sjYcs,*.) 


L 





f 


-U - ~ 




'*1 


.a. 




'J- 


(Xr-$) 


L 

(\sh (f,?U )~t Cf.*i) )Y<X ^j)J 1 f 


1: 




Of J 


C Xr ~f) 


where f is the X-coord variable, >i the Y-coord variable 
and Xr is the X-coord reference point (called the receiving 
point in formulation) 


3 . Equations and Procedures : 

In the above equations , JEL is a vector and is written (in 
rectangular coordinates) as 

JTL = [X, Y, 1] [T] 

where the area matrix T is given by 

pi X 2 til'' 

T = y i Y 2 Y 3 

(l 1 1. 

The X^, Y i relate to the node points of the triangle in 
question. If the determinant of T is 1E-G6, no calculations 



are made, and the value of zero is returned for the integral. 
A 5-point Gaussian integration technique is then applied as 
a numerical procedure. Checks are made during each step on 
the upper and lower limits to ascertain that the total strip 
area remains within the mach cone area (i.e. it may be a 
partial element configuration'). Since fl is a 1*3 array, 
the potentials returned are stored in a 1*3 array (complex). 
When possible, calculations are made in real numbers to 
computation time. A branch is provided to bypass calculation 
of integral b (above) when necessary. 


4 . Input Arguments : 

XI, X2 , X3 = ordered set of X-coord values of the 

triangle (smallest to largest) 

Z = Z-coord of triangle (Z1+Z2+Z3/3) 

A12, A13, A23 

B12 , B13, B23 

CASE 

5. Output Arguments ; None 

6. Error Returns : None 

7 . Calling Sequence : 


= slopes of sides of triangle 

= X-intercept of slopes 

= (1 or 2) relative location of the 

base of triangle 


CALL TRINTX (XI, X2, X3, A12 , A13, A23, B12, B13, CASE, Z) 


8 . Storage Required : 

9. Subroutine User : 

10 . Subroutine Required : 

1 1 . Remarks : 

A typical element 
is shown here to 
indicate what points 
are chosen for 
integration purposes 


3956 bytes (989 words) 
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APPENDIX A 


As a debugging aid, intended for the programmer only, the 
input package will accept an LABEL data section not previously 
described in the users manual. This section is called "PRINT" 
and requires 2 (two) cards. Card 2 tells the program which areas 
are to provide printout. Any or all of card 2 may be blank or 
zero. (A blank/zero inplies print off.) 

Card 1 : 

Card 2 : Format :Type 3 


1 2 3 4 5 

|P| RlIlNlT 


Format: Type 2 


MAIN 1 

INPUT, DATA, DIAPH3 1 

MASTER, ORIENT •_ ' 1 

TRINTX 2 

PSI 1 

DWASHY • 2 

ASSPOT, ASSY 1 

ASSDWN 1 

CMINV 1 

TAPES 3 

STORE, TINOUT : 1 

SIZE : 1 


This section is optional. It may appear anywhere before the 
(first) end card. (A : logical position is between TITLE and SYSTEM) 
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** NOTE ## 


Judicious use of this option is suggested, as horrendous amounts 
of print can be generated. The potential of this option to the 
programmer can be shown by exercising the full print on the very 
small test case supplied with the users manual (Ref. 2 ). 



APPENDIX B. 


PROBLEM SIZE ALTERATI ONS 

This program utilizes both fixed and variable length arrays. 

If a change in program size (data) is required, only the fixed 
length arrays need be changed. Table B.l will be the basis for 
further discussion. 

If data storage must be altered, the user must first determine 
his upper limits and then alter the arrays according to the following 
definitions. - • 


1 . XYZ (i, j, k) (grid point coordinates) 
i = coordinate value (x, y, or z) 

j = grid point number (1 — N) 
k = planform ( 1 or 2) 

2. IBLN (i, j, k) (element boolean) 
i = grid point number 

j = element number 
k = planform (1 or 2) 

3. XMOD (i, j, k) (mode shapes) 
i = mode shape for grid point 
j = mode shape set 

k = planform (1 or 2) 
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4 . LND (i, j) (leading edge grid point #) 

i = grid point number (increasing order) 

j = planform (1 or 2) 

5. MLINE (i, j) (mach line grid point #) 
i = grid point number (increasing order) 
j = planform (1 or 2) 

6. TITLE (i, j) (work storage) 
i = 4-literal characters 

j = card number 

7. INDWN (i, j) (input downwash - calculated) 

i = related grid point value (planform 1 and 2) 
j = related mode shape set 

Note: planform 1 data is stored in "i" col, position 1 through 

NP1; planform 2 data is stored in "i" col. 

( NP1 + 1) through NP = (NP1 + NP2 ) . 

Once the new problem size is determined* program alteration takes 
place in the following manner: 

Step 1. Subroutine SIZE- alter data stament (card no. 00400050) 
as required. 

Step 2. Main program - change cards' as required. 

dimension of WORC (card no. 00100010), 
dimension of WORK (card no. 00100040), 

set new maximums on cards no. 00100500 through 00100580 
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Step 3. Check Table B.l for array/subroutine and alter dimension 
statements as required. 

Step 4 . Exercise new problem data (no run card) with PRINT option 
on for sub SIZE. Compare with hand calculations. 




Delivery 

Size 

• 

1 ■ ■ . i 

Subroutine 

Array 

• 

Type 

' • 

Input 

Data 

DIAPH 3 

Master 

Output 

DISPLA 

• 

• - 

XYZ 

R 

' 

(3,80,2) 

Y 

Y 

Y 

Y 


. 

IBLN 

I 

(3,100,2) 

■ 

Y 

Y 

Y 

• ■ 

Y 


- 

XMOD 

R 

' 

( 50 ,io, 2 ) 

■ 

Y 

Y 

- 

Y 

. 

Y 

- 

LND 

1 

(20,2) 

Y 

Y 

■ 

Y 

' 

Y 

■ 

. 


MLINE 

I 

(20,2) 

. 

Y 

• 

Y 

Y 

- 

- 

TITLE 

R 

(20,20) 

Y 

Y 


■ 

- 

: - ' 

- 

- 

INDWN 

C 

(100,10) 

- 

■ - 

■ 

Y 


. Y 1 

i 


Table B.l 

Legend : 

Y = array present in subroutine 

R ~ Real 

I = Integer 

C = Complex 
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APPENDIX C 


SAMPLE PROBLEM/DATA/RESULTS 


The purpose of this section is to demonstrate graphically 
how the user describes a planform, the actual data required, and 
the results generated by the computer. Figure C.l shows a simple 
DELTA-WING planform consisting of ^-gridpoints and 2-elements. The 
broken line elements and gridpoints are generated interally, but 
the user should always have an idea of what this "DIAPHRAGM" region 
looks like before data preparation begins. The angle & «. 2,1,6) 
is computed thru SINE 0=1/ (MACHNBR). Gridpoints are indicated 
by numbers. Elements are the circled numbers. Note that the 
coordinate system is a normal right-handed system if this page is 
turned 90-degrees. The CHORD of this planform lies on the X-axis, 
with gridpoint '1 at the origin. 

Should 'questions arise during the study of this example, 
re-read the LABEL section covering the question. 


v y 



Figure C.l Sample Planform Showing Basic Structure 
and Program Generated Elements 
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COL 1 COL 2 

ROW REAL IMAG PEAL IMAG f 

*«*»<.*»**»»« ««*•**** ************** **»<*****»**»**♦***** ************** 

1 0.11337*6 00 0. 291 563E 00 0.152803E 01 -0.382149E 00 

2 0.38C8KE-01 0. 106393E 00 0.555082E 00 -0.113576E 00 



MACH NOR.. 3.U50E 01 
FrEOJENC*.. 0.4000E 00 





COL 1 COL 2 

ROM REAL 1MAG PEAL IMAG 

********* *** ************ ************************************** ****** 

1 0. 1 1 337AE 00 0. 291E63E 00 O. 1520O3E 01 -0.382199F 00 

2 0.3BC81AE-01 0. 106393E 00 0.955083E 00 -0.113576E 00 
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